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Abstract 

Continuum field control model of a street canyon is considered. The six separate 
monocriterial optimal control problems consist of minimization of functionals of the 
total travel time, of global emissions of pollutants, and of global concentrations of pol- 
lutants, both in the studied street canyon, and in its two nearest neighbour substitute 
canyons, respectively. The six optimization problems for the functionals are solved 
numerically. General traffic control issues are inferred. The discretization method, 
comparison with experiment, mathematical issues, and programming issues, are pre- 
sented. 

1 Introduction. 

The present article is a continuation of [1]. The notation of [1] will be used. Now, we will 
solve numerically six separate monocriterial optimization problems 01-06. We assumed 
the data from real street canyon of the Krasihski Avenue in Cracow [2] . We solve the set of 
nonlinear partial differential equations E1-E8 with given boundary B0-B8, and the initial 
conditions C0-C8, and sources D0-D2, by the finite difference method, using the C language 
programme written by the author. We solve this set in the cuboid Q, starting from initial 
conditions, and we iterate it over the time period using the direct finite difference method 
taking into account the boundary B0-B8 conditions, and the sources D0-D2, and initial 
conditions C0-C8, at each time step. The functionals in F1-F6 are iterated with the same 
steps as the equations E1-E8 are iterated. The first time derivative is approximated by 
the first differential quotient using forward two-point first difference in the direction of time 
coordinate, whereas the spatial first derivatives are approximated by the first differential 
quotients using central three-point first differences in the directions of spatial coordinates. 
In general, the numerical results are in good agreement with the measured data from [2]. 
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2 General Traffic Management Inferences. 



On the basis of the performed simulations, we deduce the traffic management inferences. 
In general, some vehicular traffic and hydrodynamical parameters influence the solutions of 
optimization problems 01-06, but not all of them. 

Rl. The direction of velocity of air mixture is important. The optimal pollutant concen- 
trations for both the single canyon and the canyon with its two substitute nearest neighbour 
canyons, 03, 06, are the lowest ones if the velocity components of the boundary and the ini- 
tial value problems are equal to zero; further, they are grater if the velocity has only nonzero 
vertical component v z ; then, they are greater if velocity has only nonzero x-component v x ; 
further, they are again greater if the velocity has two nonzero components (0,v y ,v z ); next, 
they are again greater if the velocity has three nonzero components (v x ,v y ,v z ); finally, they 
are the greatest if the velocity has only nonzero y-component v y . We infer that in case 
R1.6, the optimal green time g%c is different from other cases. It means that velocity of the 
mixture influences optimal vectors of control. 

R2. The optimal values for cases Ol, 04, for cases 02, 05, and for cases 03, 06, 
decrease with "uniformization" of the vehicles, when we pass from nonuniform vehicles 
UNIFORM = to uniform ones UNIFORM = 1. It is a result of reduction of the number of 
vehicles moving in the canyon. For UNIFORM = 1 the values of maximum free flow speed, 
jam, saturation, threshold, green, and red densities reach the minimal values [3]. 

R3. The long vehicular queues decrease the total travel times Ol, 04, and they increase 
both the optimal emissions 02, 05, and optimal concentrations of pollutants 03, 06 [3]. 
The decrement of total travel times Ol, 04, with long vehicular queues is a result of 
clustering of the vehicles. 

R4. The constant of temperature scale TH does not differentiate the values of optimal 
concentrations of pollutants 03, 06 in the temperature range near standard temperature 
and pressure STP conditions. However, it diminishes them even hundredfold for very high 
temperatures . 

R5. The functional form of initial and boundary conditions affects the optima. If 
they are constant, then optimal concentrations 03, 06, are two times higher than in 
the case when they are changing exponentially in space in three dimensions. The index 
FORM = CONSTANT corresponds to constant initial and boundary conditions, whereas 
FORM = EXPONENTIAL corresponds to the spatially exponentially changing initial and 
boundary conditions. For the latter case, we chose the following functional form of ini- 
tial conditions: the temperature is equal to T (x,y, z,t) = exp(— ^ — | — ^)TH, the only 
non-zero coordinates of velocities of mixture are the x-coordinates and they are equal to 
v xfl (x, y, z, t) = exp(-f -\-- c )VX, the density is equal to p (x, y, z, t) = exp(-f -\-- c )pstp, 
the concentrations are equal to Ci t o(x,y, z,t) = exp( — - — | — -)cj ; gTp, and the pressure is 
equal to p (x,y, z,t) = exp(— | — | — |)pstp- The boundary conditions B0-B8 are similar. 

R6. The presence of vehicles on both the left and right lanes is important. The optimal 
total travel times and emissions are halved in absence of vehicles on the left or right lanes 
when compared to situation when they circulate on both the left and right lanes [3]. 
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R7. The values of saturation, arrival, or jam vehicular density, and of vehicular free flow 
velocities also affect the optima 01-06 [3]. 

R8. The assumption of energy conservation equation, of the thermo diffusion effect, 
of the chemical potential and of the Grand Canonical ensemble, as well as of influence of 
gravity on intrinsic energy and on the chemical potential, drastically changes the optimal 
concentrations 03, 06, towards the measured ones [3, 4, 5]. 

R9. The value of time of simulation and of discretization in time considerably affects 
the optima (compare [3]). The values of optimal solutions 01-06, increase from tenfold 
to hundredfold. Also the optimal vectors (5-tuples) of control for 01-06, change their 
values. It is a result of cumulative effect of length of the period of simulation T s on integral 
functionals F1-F6. 

We use the following notation: yL = 0/yR = means that the lengths of all queues on 
the left/right lanes at the beginning were equal to zero, VX = is meant for zero initial 
and boundary velocity of mixture, while VX = 1 is put for the velocity of l[m ■ s -1 ]. The 
same holds for VY and VZ. If VX ^ 0, then the left lanes are leeward and the right ones 
windward. If there are no vehicles on the left or right lanes, then LON = or RON = 0, 
respectively. UNIFORM = stands for non- homogeneous (different) values of maximum 
free flow speed w vtJ , jam k vt ^ m , saturation fc„ t)Bat , threshold ^threshold, green A^green, and 
red fc„t ; RED, vehicular densities for VT = 4 types of vehicles: passenger cars, 8-ton, 12-ton, 
and 16-ton trucks, where threshold vehicular densities ^threshold are the numerical values for 
which we assume in numerical simulations that there is a vehicle on the given lane. We also 
assume that these values are the same for both the left and right lanes. 8 x ,8 y ,8 z ,8 t , are the 
steps in x, y, z, t, directions, respectively. 8^ , 8 gi , 8c 2 , 8 g2 , 5p, are the steps in C±, gi, C 2 , <?2, F, 
directions, respectively. In the boundary and initial conditions BOa, BOb, BOc, CO, we 
assume that the temperature is equal to a given constant TH. 



3 Discretization scheme. 

In order to solve numerically the set of nonlinear partial differential equations E1-E8, one 
uses the apparatus of finite differences. Firstly, we discretize the points of domain E in the 
standard way: 

(x,y,z,t) ~ (xi,yj,Zk,ti),i = 0,---,N x ,j = 0, • • ■ , N y , k = 0, • • ■ , N z , I = 0, • • • , N t , (3.1) 

i — 1 

Xi = x + — — [x Nx - xo), (3.2) 

j - 1 

Vi = yo + -^— (vNy -yo), (3.3) 

k — 1 

Zk = z + -jj—(z Nz - zq), (3.4) 



3 



ti = t + l —^(t Nt -t ). (3.5) 

Secondly, the value of given function / at point (x, y, z, t) of its domain S is approximated 
by: 

f(x, y, z, t) ~ f(x h y u z k , *,) = / iJ>fcjJ . (3.6) 
Thirdly, the first finite differences in spatial and temporal directions: 

^xfi,j,k,l = ~(fi+l,j,k,l — fi-l,j,k,l), (3-7) 
\fi,j,k,l = 2 (Aj"+l,fc,i ~ fi,j-l,k,l), (3.8) 

&zfi,j,k,i = ^(fi,j,k+i,i - fi,j,k-i,i), (3.9) 

^tfi,j,k,l = fi,j,k,l+l ~ fi,j,k,h (3.10) 

allow to approximate the first partial derivatives: 

(3.11) 
(3.12) 

(3.13) 
(3.14) 

Since the vector and scalar partial differential equations E1-E8 are of the following form: 



df(x,y,z,t) 


^ ^xfi,j,k,l 


dx 


Ax 


df(x,y,z,t) 


^ Ayfi : j t k,l 


dy 


Ay 


df(x,y,z,t) 


_ A z fij t k,i 


dz 


~ Az 


df(x,y,z,t) 


^ ^tfi,j,k,l 


dt 


At ' 



(3.15) 



hence we approximate them as follows: 



A fa A fa A A fa 

tJ i,j,k,l _ rpa/^ + fa XmJi,j,k,l £m Uijjj \ (r>-\f}\ 

~A^~* W'Vj'W'Jw AXm ' Ax m Ax n [6 b) 



Hence, we can rewrite it in following forward scheme: 

A f a A A f a 

ca fa I A -i ™/„ j. fa ' x mJi,j,k,l x m x n Jij,k,l ffl \ f Q 1 7"\ 

Ji,j,k,i+i ~ Ji,j,k,i + ^ l ■ r K X i,yj, Z k,H,Ji,j,k,n a » a a iJi,j,k,l)- \P- L( ) 

We give now the examples of dicretized equations. The discretization of continuity equa- 
tion E2 reads: 

Pi,j,k,l+l = Pi,j,k,l ~ J^ V i,j,k,l(pi+l,j,k,l ~ Pi-l,j,k,l) (3-18) 
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At y ( ^ 
~2Ay Vid ' k ' l ^ Pi,j+1,k ' 1 ~ Pi j~ 1M > 

~2Az Vi ' j ' k ' ' j ' k+1 ' 1 ~ Pi ' j ' k ~ 1 ' 1 ' 

+At Si,j,k,l 

Po 



At 



2Ax 
At 

'2Ay 
At 



Pi,j,k,i(Pi+i,j,k,i - Pi-i,j,k,i) 
Pi,j,k,i(pi,j+i,k,i — Pi,j-i,k,i) 



■) ^ Pi,j,k,l{Pi,j,k+l,l Pi,j,k-l,l)- 

For the x-component of Navier-Stokes equation El we have: 



'i,j,k,l+l 

At 



= V, 



i,j,k,l 



At 
2Ax 



v i,j,k,i( v i+i,j,k,i v i-i,j,k,i) 



2Ay 
At 



v i,j,k,l( v i,j+l,k,l 



j-l,k,l) 



' 2Az Vi ' j ' k ' ' j ' k+1 ' 1 Vi '^ k ~ 1 ' lJ 



-At- 



Tq v x 



S u iJ,k,l 



PO Pi,j,k,l 



s 



i,j,k,l 



At 



T 2 



2Ax 



Tq I ) , . I 

2 [Pi+l,j,k,l - Pi-l,j,k,l) 



a Pi,3,k,l 

At f]T s 1 x _ 
" (Ax) 2 p a 2 P,,.u [Vi+1 ' jAl Vi - 1J ' k ' l) 

At 7]T S 1 , x x x , 

(Ay) 2 po& 2 Aj,M [Vi ' j+1 ' k ' 1 ~ W + Vi ^ l) 

At ^ 1 ( v * _2v x + v x ) 



(Ax) 2VS ' 3 7 p a 2 p iJtkjl 
At j^ + V\V T s 1 



i v i+i,j,k,i 2v i,j,k,i 



,j,k,l) 



-1,3 



(v? 



(3.19) 



V i+l,j-l,k,l V i-l,j+l,k,l + V i-l,j-l,k,l) 



4AxAy V ^3Vo& 2 PM,M V 

At V^vTs 1 / x x x x \ 

4AzAl + 3 j p c 2 A,i,fe,/ {Vi+1 ' j ' k+1 ' 1 ~ Vi + 1 ^ k - 1 ' 1 ~ v ^,k+i,i + v i-i,j,k-i,i) 



+At 



9oT$ g * 
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Similarly for the y-component: 

u i,j,k,i+i — u i,j,k,i 2Ax ' j,k ' +1 ' j ' k ' 1 u i-i,j,k,i) 

~2ty V ij,k,l( V h+l,k,l ~ V i,j-l,k,l) 

~2^^ V i,j,k,l( V i,j,k+l,l ~ V i,j,k-l,l) 
a.Ts V i,j,k,l q 

-At—- bij tk ,i 

Po Pi,j,k,l 

At T| 
--^jBT — (T ijj+hk) i - Tij-^kj) 

3,k,l i \ 
BTq — [Pi,j+l,k,l - Pi,j-l,k,l) 



2Ay b 2 p iJikil 

At 1 ( v v _ 2v y + v y 

(Ax) 2 p a 2 p iJiktl 1 Vi+ W ^ + 

At 7]T S 1 



(Ay) 2 p b 2 p iJjktl 
At r,T s 1 „ 



i,j+l,k,l ^ v i,j,k,l + v i,j-l,k,l) 



V i,j,k+l,l 2v tj,k,l + V i,j,k-l,l) 



(3.20) 



(Az) 2 p c 2 p iJA i 

At (3 \ 7] ) r]Ts 1 (v x -v x -v x +v x ) 

AAxAy ^ + 3 j p b 2 p i:jAl V *+U-W V i-U+W + v i-U-hk,i) 

At i] r]T s 1 



^ L tc i 'L\ ']±JL 1 i y y 



i-v 



i-ij+i 



,k,l + v i-l,j-l,k,l) 



At (3 \ V ) vTs 1 (v x -v x -v x +v x ) 

AAyAz 3 p b 2 p ijk i ^ { >i +1 ' k+1 > 1 v iJ+hk-i,i v i,j-i,k+i,i + v i,j-i,k-i,i) 

+At^flg y . 
b 

Finally, for the z-component: 



v lj,k,i+i - v tj,k,i 2A.x V *' j ' k ' l ^ V * +1 ' j ' k ' 1 v ^~ 1 '3' k ^ 



(3.21) 



At 



2Ay 
At 



v lj,k,i( v i,j+i,k,i v i,j-i,k,i) 



2Az Vi ' j ' k ' ' j ' k+1 ' 1 Ui >j< k - 1 ' lJ 

i,j,k,l 



-v 



_£ dZ s_jjM s 



PO Pi,j,k,l 

At 



2Az 
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At Tij kj 

oaT °72"7 ±J - L {Pi,j,k+i,i - Pi,j,k-i,i) 

ZLAZ C Pi,j,k,l 

At rjT s 1 



(Ax) 2 p a 2 Pi, j: k,i 
At r]T s 1 ( 

(Ay) 2 p b 2 p iJA i 
At t?T 5 1 ( 

(A2:) 2 p c 2 p iJjfc> , 



( v i+l,j,k,l ~ ^ V i,j,k,l + V i-l,j,k,l) 
\ V i,j+l,k,l ~ ^ V i,j,k,l + V i,j-l,k,l) 
, V i,j,k+l,l - ^ V i,j,k,l + V i,j,k-l,l) 



At _^ + V^vTs 1 



4AxA2; 3 p c 2 pij, k ,i 
At , 77 1 



AAyAz 3 p b 2 p iJ)kt i 



^ v i+i,j,k+i,i v i+i,j,k-i,i v i-i,j,k+i,i + v i-i,j,k-i,i) 



, v i,j+l,k+l,l V i,j,k-l,l v i,j,k+l,l + V i,j-l,k-l,l) 



At 1] 7]Ts l/ 2 z 2 

(s + ~ (^ij.fe+l,/ _ 2t MVM + V i,j,k-l,l) 



(Az) 2 
+At 9 -^- 
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4 Method of comparison. 

We briefly describe the method of comparison of the numerical results with the experimental 
ones. The comparison of simulations with the measured data, if authorized, is aimed at 
showing the correctness of the model (conservation of the rank of simulated parameters), 
and of presenting the correct direction of deviations (the deviations always tend to secure 
the direction). We transform both measurements and optima Jq 03, in order to compare 
them. The measured pollutant concentrations are multiplied by the volume of the canyon, 
and we obtain pollutant masses (we assume that the measured concentrations are the same 
in the entire canyon and in time). The pollutant masses are added up, assuming HC mass 
to be equal to zero (HC were not measured in Ref. [2]). The calculated optima Jq 03 
are divided by time of simulation T5, giving approximation of pollutant masses in every 
moment of simulation (we assume that Jq 03 are homogeneous in time, in order to make 
possible the comparison with measurements). Thence, we have two comparable magnitudes. 
Also, only insufficient vehicular flow data are given in Ref. [2]. Neither the jam, saturation 
vehicular densities, nor the maximum vehicular speeds are present. However, the initial 
and boundary conditions of traffic flow B5-B8, C5-C8 are not constant functions because 
they are parametrized by measured traffic parameters. The data are not instantaneous, but 
they are measured in intervals of hours. Moreover, the measured data are not accompanied 
by measurements of wind velocity, temperature, density. If the measurements had been 
performed instantly using more specialized equipment (e. g., lidar), and if they had been 
stored in computer facility, then, it would be possible to make more serious comparisons. The 
simulations were performed for the time period T = 60 [s], whereas the measurement were 
averaged over the time period T = 1800 [s]. The integral functionals F3, F6, are cumulative 
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with respect to the time period. 



5 Computer simulations. 

From thousands of performed optimizations 01-06 we selected some representable opti- 
mizations. The C language programme of more than 14000 lines written by the author was 
being run on workstations of Hewlett-Packard, Sun, Silicon Graphics, under UNIX operating 
system. It was being run also on Pentium personal computers under Linux operating sys- 
tem. We used C language UNIX/Linux compilers cc. The time of simulations varied from 
a couple of hours to several days or even weeks, depending on the discretization parameters. 
The discretization steps were tested in order to obtain finite solutions to the optimization 
problems in reasonable time of simulation. They were also studied from the point of view of 
sensitiveness of optima appearing on them. 

6 Mathematical questions. 

The problem of uniqueness of solutions of optimization problems 01-06 is complex. If 
there are no vehicles allowed in the canyon, then all the six functionals 01-06 are constant 
functions of parameters, the optima exist, and the number of optimal solutions is continuum. 
The values of functionals are then zero for total travel time and emissions, and are positive 
for concentrations. These results are analytical. For other cases there are no such analytical 
results for both the uniqueness and existence known to the author. The numerical solutions 
of optimization problems are approximately global within the error induced by discretization 
of the physical domain E and of the domain of control parameters f/ adm . The method of 
optimization was a full search. 

7 Conclusions. 

The proecological traffic control idea and advanced model of the street canyon have been 
developed. It has been found that the proposed model represents the main features of 
complex air pollution phenomena. 
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